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ABSTRACT 

The Finite Element Analysis Program (FEAP) was modified and 
integrated with the Naval Postgraduate School version of the 
Computer Analysis Language (CAL-NPS) . This enables the 
solution of linear and non-linear, two and three dimensional 
heat conduction problems in an interactive mode. The usual 
types of boundary conditions, including radiation, may be 
specified. The heat conduction group includes prompts for 
user supplied data. Several existing CAL-NPS commands were 
improved and a ’'HELP” facility was added. Commands were 
added for visual display of the finite element mesh at 
graphics terminals. The User Guide for this expanded 
version of CAL-NPS is provided. 
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I. IHTRODacnOH 



A. GENERAL DESCRIPTION 

Since the implementation of the Finite Element Analysis 
Program (hereafter ref ere d to as FEAP) at the Naval 
Postgraduate School, it has had only limited use. It is now 
available to general NPS users in an interactive form 
through the IBH 3033 VM/MVS time sharing system. The data 
management and user interactive command structures were 
established within the existing interactive program. 
Computer Analysis Language (CAL). This integration of 

systems provides the ability to solve linear and nonlinear, 
steady and unsteady, two and three limensional heat conduc- 
tion problems involving temperature dependent thermophysical 
properties and complicated radiation/convection boundary 
con ditions. 

Additional capability was provided to CAL system users 
both in changes to existing subroutine groups and in the 
addition of a new graphics group with its attendant command 
structure. The graphics functions enable the user to plot 
two- and three- dimensional structural and heat transfer 
meshes. Another new facility is the HELP operation which 
allows a user experiencing trouble with a particular opera- 
tion to interactively obtain assistance. 

B. HISTORICAL BACKGROUND 

The original CAL program was developed by Professor E. 
L. Wilson of the University of California in 1977 [Ref. 1 ]- 
It was later adapted and modified for use at the Naval 
Postgraduate School by LCDR L. B. Elliot [Ref. 2]. The FEAP 
program -was written by Professor R. L. Taylor of the 
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University of California [Ref. 3] in 1977. Implementation 
of FEAP at Naval Postgraduate School was done by LT J. H. 
Bettencourt [Ref. 4]. 

C. OBJECTIVES 

The objectives of the author's work have been to: 

1. integrate the data management system of FEAP with CAL 
to create an interactive conduction heat transfer 
problem solving system; 

2. modify existing operations to extend their useful- 
ness; 

3. add operations to extend the capabilities of the 
program to include graphics; 

4. create a HELP facility; 

5. create a USER'S MANUAL to facilitate use of this 
program. 
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II. ORGANIZATIOR OP HEAP rRANSPEH GROUP 



This chapter provides a general overview of the organi- 
zation of the Heat Transfer Group of CAL. It is intended to 
provide sufficient information to permit users to operate 
the CAL heat transfer package. 

The execution of the program is flexible and controlled 
by user selection of operations in a logical sequence from 
the commands that are available. There are two broad 
categories of operations, data input and problem solution. 

A. DATA INPUT 

1 • In i tializatio n 

The heat transfer group performs matrix creation and 
manipulation automatically. As the problem progresses, 
arrays are created, altered and deleted under program 
control. Through the HTXFR operation the user provides 
sufficient information to establish the initial arrays for 
data input and problem solution. The number of nodes, 
number of elements, number of material sets, spatial dimen- 
sion, number of degrees of freedom per node and the maximum 
number of nodes per element are required. For heat transfer 
problems the number of degrees of freedom per node is always 
one. The option to assign a higher number is available 
because the equation solvers in this program are applicable 
to other fields of which future work may make use. 

2 . Nodal C oordinate s 

Nodal coordinates are input via the COORD operation. 
This operation has built in node generation capability. By 
specifying an initial point and a node generation vector. 
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the user may easily input large meshes. Coordinate system 
conversion is also available. Coordinates may be input in 
the Cartesian system, the cylindrical system with any one of 
the three axes longitudinal, the spherical system or any 
combination of the above systems. All coordinates are 
converted to Cartesian coordinates for use with CAL. 

3- Element C onnect! vi t y 

The ELCON operation inputs the element connectivity 
data. Here again is a generation capability. The user may 
specify the connectivity for one element and a generation 
vector to create additional rows. 

For the two-dimensional elements the user may 
specify a 4 to 9 node isoparametric element. There is an 8 
to 21 node isoparametric element for three-dimensional 
elements. Both of these elements must follow the numbering 
convention shown in Appendix A. 

^ • Hater ial Propertie s 

The required amount of material property information 
varies from problem to problem. The PROP operation prompts 
the user for the information required to solve the problem 
at hand. At a minimum the material's conductivity (k) , 
specific heat (c) , specific mass (rho) ,heat generation per 
unit volume (q''') and the geometry type (plane or axisyme- 
tric) must be specified in a consistent system of units. 
Appendix A includes examples of consistent systems. 

The user also inputs the number of Gaussian points 
per direction for quadrature and codes for temperature 
dependent properties and boundary oonditions. The codes 
indicate other information required. Temperature dependent 
properties are input as tables and Linear interpolation is 
used to determine the property value at a given temperature. 
Boundary conditions are specified as shown in Appendix A. 
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If an exterior boundary line or surface condition is not 
specified, it is assumed to be insulated. 

5. Constant Temperatu r e No de s 

For problems requiring certain nodes to be at 
constant temperature, the CTEMP operation is available. 
This operation may also be used during the solution stage of 
the problem to provide step changes at previously specified 
constant temperature nodes. Because this operation gener- 
ates an array used in profiling the solution equations, the 
user may not change the node numbers that were established 
as constant constant temperature nodes after execution of 
the PROF operation. The temperatures of these nodes, 
however, may be changed. 

6 . E quation Profile 

The PROF operation establishes the equation profile 
for problem solution. Prior to the execution of this opera- 
tion, any data input may be changed by specifying the appro- 
priate operation and re-entering the data. After its execu- 
tion the user may not change the nodes designated as 
constant temperature nodes to temperature varying nodes or 
vice versa. The user may change the value of the constant 
tern pe natures. 

B. SOLOTIOH 

The matrix formulation of the heat transfer problem as 
discussed in Reference 4 is: 

tK]CT} + [C]{T} + {F} = CO} (1) 

where (K) represents the conductivity matrix, (C) represents 
the heat capacity matrix, (F) represents the flux vector and 
(T) represents the temperature vector. The derivative of 
{T} with respect to time is {T} . The flux vector includes 
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heat generated per unit volume and boundary fluxes as speci- 
fied for given boundary surfaces. This is a fully general- 
ized formulation, including non-linearities, since the 
matrices [K] and [C] and the vector {F} can be temperature 
dep en dent. 

1 • F orming Cond u ctivi t y Matrix 

For heat transfer problems, which involve only one 
degree of freedom per node, the conductivity matrix, [K], 

will always be symmetric. The capability for generating 

unsymmetric matrices was provided, but will only be appli- 
cable when additional types of problems are programmed into 
the CAL system. The command for the unsymmetric conduc- 
tivity matrix formulation is OSYMC. 

2« Forming Heat Capac ity Matr ix 

The heat capacity matrix, [C], used in time depen- 
dent problems, can be generated with either CCAP or LCAP. 
To form a consistent capacitance approximation use the oper- 
ation CCAP. A lumped capacitance approximation is formed 
using the LCAP operation. 

3. Forming Flux Vecto r 

To complete the problem formulation, the flux 
vector, {F} , must be generated. The FORM operation forms 
the flux vector taking into account the internal heat gener- 
ation and boundary surface fluxes as indicated in the PROP 
operation. 

4 . Eg n ation Solvin g 

Once the time independent problem is formulated, the 
temperature vector, {T} , is calculated by the CALC opera- 
tion. Time dependent problems do not use this operation, 
but rather the ordinary differential equation solver. 
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5. First O rder Ordinary Diff erential Equat ion Solver 



The first order ordinary differential equation 
solver is accessed with the ODE operator. It employs the 
Zienkiewicz two- and three- level schemes [Bef. 4, 5]« 

In addition to tha time step size change using the 
DTIM operation, an optional automatic time step adjustment 
is incorporated in the ODE operation. The norm of the 
difference between temperature vectors at two consecutive 
times is computed at each step. If the norm is less than a 
user specified maximum temperature difference, the time step 
is doubled before going to the next step. If the norm is 
greater than a user supplied minimum temperature difference, 
the time step is halved and calculation for that time step 
is repeated until the norm is acceptable. If the tempera- 
ture differences are specified as zero, no time step adjust- 
ment will be performed. 

The user specifies one of three functions which are 
performed by the ODE operation. A second operation name, 
which must be separated by a comma, follows ODE. The 
options are INIT, LINE or aOAD. 

The operation ODE, INIT is used to input the integra- 
tion constants theta, beta and gamma [Ref. 5], the maximum 
and minimum temperature differences for the automatic time 
step adjustment and the initial temperature vector. No time 
integration is performed by this instruction- 

ODE,LINE performs the two point scheme and the 
current temperature vector is substituted by the newly 
calculated temperature vector. 

The ODE, QUAD operation is similar to the ODE, LINE 
operation but uses the three point scheme. 

6- Pr inting No^l Temp eratures 

Once the temperature vector has been updated, the 
PTEMP operation prints the temperatures in node number 
order. 
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III. ORGANIZATION OP GRAPHICS GROOP 



The graphics group is capable of displaying meshes 
either on the IBM 3277 dual screen terminal system or any 
PLOT-10 compatible terminal. It is initiated through the 
use of the GRAPH operation by which the user specifies the 
type of graphics terminal in use. 

A. TITLE 

The TITLE operation is used to title the mesh being 
displayed. The user may input up to three lines of fifteen 
characters for the title. This operation must be specified 
immediately prior to the displaying operation. 

This operation calls the USRIN subroutine which reads 
from the terminal three lines of characters. A flag (IFLAG) 
is set to indicate a title is to be printed. 

B. HEAT TRANSFER MESH 

The PLHX operation locates the heat transfer coordinate 
and connectivity arrays and displays the mesh. The viewing 
area is optimized so the longest dimension is full screen. 
The maximum and minimum coordinate values for each direction 
are displayed. This operation displays two- and three- 
dimensional meshes. If the mesh is three-dimensional, the 
user must specify the viewing plane. 

This operation calls either the FPPLOT (2-D) or FP3PLT 
(3-D) subroutine. These subroutines scan the coordinate 
arrays for maximum and minimum values and initialize the 
graphics screens. They both use subroutine BOX to set the 
virtual window, set the screen window and draw a box around 
the plotting area. The standard element connectivity is 
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stored in an array which is used in conjunction with the 
user’s element connectivity matrix to draw the elements line 
by line. The center node is plotted using a symbol. 
The user’s title is plotted (if IFLAG is not 0) by subrou- 
tine OSRTIT, which then resets IFLAG. The maximum and 
minimum values of the vertical and horizontal coordinates 
are printed by SCRDAT. Subroutine SCR DAT also indicates axis 
orientation. Prior to terminating the screen graphics, 
subroutine TITLE is used to write the title box identifying 
the type of mesh being displayed. 

C. STHUCTORAL BESH 

The PLST operation displays the structural system speci- 
fied by user supplied coordinate and connectivity arrays. 
The viewing area is optimized as described for the PLHX 
operation. Likewise, the user may specify the viewing plane 
as either the X-Y, Y-Z or X-Z plane. 

This operation calls the CLPLOT subroutine. CLPLOT uses 
the same logic and subroutines as FPPLOT. The difference 
between the two subroutines is in the method of storage of 
the coordinate arrays and the connectivity matrix. 
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IV. CHANGES TO EXISTING GROUPS 



Modifications were made in several existing subroutine 
groups to make their operations more versatile. One utility 
subroutine was improved and one was added. 

The changes that were made are sensitive to prior 
versions of CAL. The same results are obtained for previ- 
ously existing operation commands. No files used with other 
editions of CAL need to be modified to operate with this 
program. 

A. UTILITY SUBROUTINES 

1 , Subroutine RCARD 

Subroutine RCARD reads and interprets the operation 
commands. An operation command has the form 
OP,M1rM2,M3,M4,M5 ,N1, N2,N3,N4 

where OP is the operation name, M1 to M5 are matrix names 
and N1 to N4 are integers. Previously, N1 to N4 had to be 
values greater than or equal to zero. The symbol/state 
logic matrix and subsequent action codes were modified to 
allow users to input negative integer values. 

This is important , for the ZERO operation which 
enables users to create matrices with a given value in the 
diagonal locations and another value in the off diagonal 
locations. Prior to this change, a user desiring to create 
such a matrix with negative values had to input the matrix 
row by row. 

2 . Sub routine FRTC MX 

Subroutine FRTCMX is a new utility subroutine which 
allows the CAL program to invoke most CP/CMS commands. 
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After the invoked command is executed, control returns to 
CAL. 

It is presently used in the SAVE and RESUME opera- 
tions to invoke the CMS command 

FILEDEF NSAVE DISK Ml SAVE (RECFM VS LRECL 7290 
BLKSIZE 7294) 

where NSAVE is a logical unit number assigned to the SAVE 
and RESUME operations and M 1 is a user input file name. 

B. GROUP 1 

1 • SA V E Operation 

The save operation creates a file on the user's 
A-disk containing all arrays in storage at the time of 
issuance. Previously the entire 100,000 word main array was 
stored, regardless of how many words were actually being 
used. It was always stored under the name FILE 02, 

preventing the user from saving more than one problem. 

The method of storing the array was changed to store 
just the number of locations actually being used. More than 
one problem may be saved because the subroutine FRTCHX was 
used to create a SAVE file with a name assigned by the user. 
If a name is not specified FILE 02 will be the name of the 
saved problem. 

2 . RESUME Operat ion 

The RESUME operation reads a saved file into memory. 
It was altered to read named files saved by the new SAVE 
operation. If a name is not specified FILE 02 will be read. 

3. LOAM Operati on 

The LOADI operation loads integer arrays. The arays 
were input row by row. The option to generate arrays was 
added. The user may specify one row and a row generation 
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vector. The number of rows specified will be automatically 
generated. This operation was moved from the static anal- 
ysis subroutine group to the general matrix command group. 

4 , PRINT O peration 

The PRINT operation prints an array in matrix 
format. Previously it could only print arrays containing 
real numbers, resulting in the erroneous printing of arrays 
containing integers. The user may now specifiy whether the 
array to be printed contains real or integer values. 

5 • BILP Operat ion 

The HELP operation was added to the general command group 
and provides the user with information on the use of all the 
available operations. It accesses a file of instructions, 
sorts through them to find the desired operation and 
displays the appropriate information on the screen, 

C. GROOP 2 

1 • NODES O peration 

The NODES operation creates the matrix of nodal 
coordinates for a structural problem. This information was 
entered in cartesian coordinates, node by node. 

The user may now opt to enter data in cartesian, 
cylindrical (any axis longitudinal) or spherical coordi- 
nates. It will be converted to the cartesian coordinate 
system used by CAL. Additionally, the user may generate new 
nodes by specifying one node and a node generation vector. 
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V. SOlPTION OF CONDOCTION HEAT TRANSFER PROBLEMS 



In order to solve a conduction heat transfer problem, 
the user must provide the solution algorithm to CAL. A 
discussion of the matrix manipulation and equation solving 
techniques can be found in Reference 4 . This chapter will 
present possible algorithms for solution of this class of 
problems. 

A. STEADY STATE PROBLEMS 

These problems take the form: 

[K]{T} ♦ {F} = {0} (2) 

1. Linea r Heat C onduc t ion Probl e ms 

This is the simplest case to consider. After input- 
ting mesh data, the user must form the conductivity matrix 
[K], (SYMC) and the load vector {F} , (FORM). Then the nodal 
temperatures must be calculated (CALC) and printed (PTEMP) . 

Consequently the sequence of solution operations for 
this type of problem would be: 

SYMC 

FORM 

CALC 

PTEMP 

2. Non-linear He^ Con duct ion Pr oblem s 

Since the conductivity matrix is time dependent on 
temperature, an iterative algorithm must be used. This 
requires a looping operation (LOOP, NEXT) around the linear 
steady state sequence. 
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The solution operations may be: 

LOOP,N1 

SYMC 

FORM 

CALC 

PTEMP 

NEXT 

where N1 is the user's guess of the numbex of iterations 
necessary to obtain equilibrium, Howeverr the program main- 
tains an internal check on the residuals. When they 
decrease below the predefined tolerance (TOL command, 
default is 10”^) , the looping operation is terminated upon 
the subsequent NEXT command, 

B. TINE DEPENDENT PROBLEHS 

These problems involve the full form of equation (1) : 
CK]{T) + tC Ht} ♦ (F) = {0} (1) 

1. Linea r Heat C onduc t ion Probl e ms 

This case requires the solution of a first order 
ordinary differential equation (ODE), Additionally a heat 
capacity matrix, [C], must be formed (CCAP or LCAP) and a 
time step provided (DTIM) , 

The differential equation solver is accessed using 
the operation 0DE,m where M1 is one of the following: 

INIT to specify initial temperature vector and 
the integration constants 

LINE to perform the two-time level algorithm 
QOAD to perform the three-time level algorithm 
The heat capacity and conductivity matrices are 
unchanged in a linear problem and must be placed outside the 
loop. The load vector is reformed every time step and the 
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These operations are 



time must be advanced with ADTIM. 
included in the loop. 

The sequence of operations to solve a linear time 
dependent problem may be: 

DTIM 

ODE, I NIT 
SYMC 

CCAP (or LCAP) 

L00P,N1 

FORM 

ODE, LINE (or QUAD) 

ADTIM 

PTEMP 

NEXT 

where N1 is the number of time steps the user wants to take. 

2- No n -linea r Heat Con duct ion Pro blems 

The heat capacity and/or conductivity matrices are 
temperature dependent in this class of problem. The temper- 
ature dependent matrix (matrices) must appear inside the 
looping operation whereas the constant property matrix would 
be excluded from the loop. 

To solve a fully non-linear problem the following 
operation sequence may be used: 

DTIM 

ODE,INIT 

L00P,N1 

SYMC 

CCAP (or LCAP) 

FORM 

ODE, LINE (or QUAD) 

ADTIM 

PTEMP 

NEXT 

where N1 is the number of time steps the user wants to take. 
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C. NOHEBICAL EXAMPLES 



1 , Hollo w Cyl i nder with Circumfe r ential Hea ting S trip s 

A hollow cylinder with a four inch outer diameter 
and a three inch inner diameter was subjected to an axial 
forced connvection condition in a wind tunnel by Professor 
P. F. Pucci of the Naval Postgraduate School. There were 
sixteen one-quarter inch wide heating strips equally spaced 
over 180 degrees of the outer surface as illustrated in 
Figure 1. 

The heating strips were maintained at a constant 
temperature of 160® F and measurements of the surface temp- 
eratures between the strips were made using teledeltos 
paper. The ambient temperature was 60® F. The cylinder was 
considered to have a density of 70 Ib/ft® and a specific 
heat of 0.6 BTU/lbm®F. Tests were made with four heat 
transfer coefficients (h). 

The finite element model took advantage of the 
cylinder’s symmetry, consisting of one half of a heating 
strip plus one half of the interval between strips. The 

element mesh, as generated by the PLHX operation, is shown in 

Figure 2. The comparison of the model data to the tele- 
deltos paper measurements was very favorable and is shown in 
Figure 3. 

2» Transient Sur fac e T empe nature s in ^ Inf in it e Pla te 

An infinite flat plate was considered as a test 

problem (Figure 4) . The plate was at an initial uniform 

temperature greater than the ambient temperature and then- 
exposed to convection conditions with a constant external 
heat transfer coefficient (h) . 

The approximate solution was obtained using a 
Heisler chart. The temperature of the outer wall was 
computed at 1 minute, 5 minutes and every 5 minutes there 
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after up to 30 minutes. The comparison between these values 
and those generated by CAL was very close and is shown in 
Figure 5. ' 
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71. CONCLUSIONS AND RECOMMENDATIONS 



The code that was integrated into CAL provides an accu- 
rate and reliable means for solving a variety of conduction 
heat transfer problems. The system is user friendly both in 
prompting for input and detecting errors. The use of the 
heat transfer group of commands is encouraged, as well as 
efforts to increase its versatility. 

While the capabilities of the program are significant, 
there is room for improvement. The present version uses 
primarily an in-core solution technique, which restricts the 
problem size to within the user's virtual machine space. 
The capacity to handle larger problems may be increased 
through the use of an out of core technique for building and 
storing the conductance and capaitance matrices, as well as 
equation solving. 

The graphics package could be expanded and improved to 
provide more information to the user. The plotting of 
isoparametric elements needs to include the capability to 
generate and plot curved lines. Another desirable capa- 
bility is to portray all three axes on the screen with 
apparent depth and allow rotation to any desired view. A 
capabilty to plot isotherms could be added to enhance the 
interpretation of the results. 
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Figure 1 . 
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SURFACE LOCATION (1/84 In) 
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Figure 4. 
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TIME (mtn) 



APPENDIX A 



OSER'S GUIDE 

This appendix provides details on the use of CAL with 
the IBM 3033 computer at NPS. The program as modified at 

the Naval Postgraduate School is subsequently referred to as 
CAL-NPS. Section A provides the details on the command 
structure. section B is a summary of commands available. 
Section C provides the job control language for executing 
the program in both the batch and interactive modes at NPS. 
Section D contains detailed specifications for each avail- 
able command. Finally, section E gives direction for 
solving larger problems with CAL-NPS. The majority of this 
appendix was originally published as Reference 1 . The 
author wishes to express- appreciation to Professor Wilson 
for permission to use this material. 

A. FORM AID RESTRICTION OF THE LANGUAGE 

CAL-NPS is an interpretive language which is designed to 
manipulate arrays and matrices, to perform standard struc- 
tural analysis operations and to perform conduction heat 
transfer analysis operations. A CAL-NPS program run 
involves the reading of the input deck once and executing 
the commands designated by the operation cards as they are 
encountered. Looping operations allow a sequence of 
commands to be executed more than once. 

The input deck is composed of operation cards and data 
cards. The data cards directly follow each operation card 
which requires data (see LOOP operation for exception to 
this) . The operation card contains the name of the opera- 
tion to be executed, names of arrays associated with the 
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Examples of the general 



operation and integer constants, 
form of this card are: 

0P,H1 ,H2,M3,M4,M5 ,111, N2,N3, N4 comments 
OP,m ,N1, N2 

op,m 

OP 

in which OP is the name of the operation to be executed. Hi 
is the name of an array and Ni is an integer. The names of 
OP or Mi are one to eight alphabetic or numeric characters 
to be selected by the user. The first character of a name 
must be alphabetic. The sequence of terms OP, Mi and Ni 
must be separated by commas. Characters following a blank 
will be printed as comments in the output from the program 
run . 

If an operation attempts to load or generate an array 
which previously existed, the program will delete the array 
before the execution of the operation. A new array need not 
be the same size of the old array which had the same name. 
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B. SOMflABY OF COHflANDS 



1 . General Commands 

indicates a significant change or addtion in CAL-NPS 



START 


- Initialize for the next problem 


STOP 


- Normal termination 


NO 


- Temporary suppression of output 


YES 


- Restores output 


LAB EL 


- Print comments 


READ 


- Change logical device for input 


WRITE 


- Change logical device for output 


TIME 


- Suppress time printout 


SAVE 


- Interrupt a problem ^ 


RESOME 


- Continue an interrupted problem ^ 


LIST 


- List arrays and- storage used 


HELP 


- Access HELP files 


2. 


General Matrix Comm ands 


LOAD 


- Load user defined real matrix 


LOAD I 


- Load user defined integer matrix 


ZERO 


- Create null or unit matrix 


PRINT 


- Matrix print operation * 


DUP 


- Matrix duplication 


ADD 


- Matrix addition 


SOB 


- Matrix subtraction 


HOLT 


- Matrix multiplication 


TRAN 


- Matrix transpose 


SCALE 


- Multiply a matrix by a scalar 


SOLVE 


- Solution of linear equations 


DOPSM 


- Form sub- matrix from large matrix 


STOSM 


- store sub-matrix in large matrix 


DDPDG 


- Form row matrix from diagonal 


STODG 


- Store row on diagonal 
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MAX 


- Evaluate row maxi mums 


NORM 


“ Evaluate matrix aorms 


IN7EL 


- Invert each term in matrix 


SQREL 

LOG 

PROD 

DELETE 


- Square root of each term in matrix 

- Natural log of each term in matrix 

- Evaluate product of all terms in a matrix 

- Delete matrix from storage 


3. 


Static Analysis Operations 


NODES 


- Input structurcil joint geometry 


BOUND 

BEAM 


- Specify boundary conditions 

- Form 3-D beam stiffness matrix 


TRUSS 


- Form 3-D truss stiffness matrix 


PLANE 

SLOPE 


- Form 3 to 8 node plane stiffness matrix 

- Form stiffness matrix from slope/deflection eg. 


FRAME 


- Form 2-D frame stiffness matrix 


LOADS 


- Form load vector 


ADDSF 

ADDK 


- Form global stiffness and mass matrices 

- Add element matrix to global matrix 


MEMFRC 

DISPL 


- Calculate element forces from joint displacements 

- Print joint displacements 


FORCE 


- Evaluate and print member forces 


4. 


Dynamic Analv sis Operations 


FUNG 


- Generate equal interval time function 


STEP 

EIGEN 

DYNAM 

PLOT 


- Integrate dynamic equilibrium equations 

- Evaluate mode shapes and frequencies 

- Evaluate uncoupled equatioas of motion by 
mode superposition method 

- Line printer plot of joint time history 


5. 


Heat Transfer Opera tions 


HTXFR 


- Initiate heat transfer problem 


COORD 


- Input nodal coordinates 


ELCON 


- Input element connectivity matrix 
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PROP - Input material property data * 

CTEHP - Input constant temperature node data ^ 

PROF - Establish profile of equations * 

SYHC - Create symmetric conductance matrix * 

OSYMC - Create unsymmetric conductance matrix * 

LCAP - Create lumped capacitance matrix * 

CCAP - Create consistent capacitance matrix ^ 

FORM - Create flux vector 

CALC - Solve time independent systems of equations * 

ODE - Solve time dependent systems of equations ^ 

PTEMP - Print nodal temperatures * 

TOL - Set solution convergence tolerance ^ 

CONV - Perform temperature convergence test 

DTIM - Set time step increment ^ 

ADTIM - Advance time by one time step ^ 

PROMPT - Suppress/restore prompts * 

6 . Graphic s Operation s 

GRAPH - Initiate graphics * 

TITLE - Label mesh plot * 

PLHX - Plot heat transfer mesh * 

PLST - Plot structural mesh ^ 

7 . Loo£ Operations 

LOOP - Start of loop 

NEXT - End of loop 

SKIP - Conditional skip of operations within loop 

8. Na mes A vailable for Use r Subr o utines 

USER A 
USERB 
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C. EXECOTIOH OF CAL-NPS 



For the time sharing (CMS) system at Naval Postgraduate 
school, do the following: 

_(Ose the standard LOSON procedue) 

_lin)c 0040P 191 199 
ENTER PASSWORD: 

_XXXXX (ESAN) 

R; T=0.01/0.01 11:09:55 

_access 199 C (Note: You must access the ”C” 

C (199) R/0 disk) 

R; T=0.01/0.01 11:10:04 
_cal 

ENTER TERMINAL CODE; 

1 = PLOT- 10 Compatible Terminal (GRAPHICS) 

2 = IBM 3277 DUAL SCREEN (GRAPHICS) 

3 = Any Alpha Numeric Terminal (NO GRAPHICS) 

_1 (or 2 or 3, as appropriate) 

(The computer responds with several lines of procedure) 
EXECUTION BEGINS 

0.0 SECONDS 

(You are now under the control of CAL-NPS) 

start 

«»*START 

! 

(Your own CAL-NPS program is inserted here) 

stop 

♦‘S'STOP 

R; T=0.01/0.01 11:12:45 
_log (Terminates session) 
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D. CAL-NPS COHHAND SPECIFICATIONS 



1 • General Matri x Oper atio ns 

CAL-NPS has most of the standard matrix operations 
plus some special array operations which are useful in engi- 
neering analysis. The following is a list of approximately 
32 operations which are used for control and general matrix 
manipulation. 

A indicates the formation of a new matrix. A 

matrix previously defined with the same name will be 
deleted. A indicates modification of an existing 

matrix. 

Note ; Whenever the expression "card” is used it is 
meant to also stand for "instruction" in interactive mode. 



START 

This operation eliminates all arrays which were previously 
loaded or generated. 



STOP 

This operation causes normal termination of a CAL-NPS 
program. 



NO 

YES 

These operations are used to selectively suppress output 
from CAL-NPS. The NO operation suppresses all printing, 
except diagnostics, until the operation YES is encountered. 
Therefore, in subsequent runs of the same CAL-NPS program, 
output which was previously correct need not be reprinted if 
these cards are inserted in the data deck. 



LABEL,N1 

This operation wj.ll read and print N1 comment cards which 
follow the operation card. Column 1 of each card will be 
interpreted as a standard carriage control symbol (i.e. 0 

for double space and 1 for skip to the top of the next 
page) . 



39 



READ,N1 

THIS OPERATION IS VALID ONLY WITH THE CP/CHS TIME SHARING 
SYSTEM. 

This operation permits the selection of a user's file. or the 
terminal as the input file device. The default is the 
terminal. If N1 is 4 , subsequent commands will be read from 
FILE FT04F001 on the user's A-disIc. If N1 is 5, the 
terminal will be restored as the input file device. All 
disk files prepared for use with this command should end 
with either STOP or READ, 5. This command will not be 
executed on the MVS (batch) system. 



WRITE, N1 



THIS OPERATION IS VALID ONLY WITH THE CP/CMS TIME SHARING 
SY?TEM, 



This operation permits the 
A-disk or the terminal as 
default is the terminal 
printed at the terminal rag 
selected. If N1 is 8, subs 
FILE FT08F001 on the user 
terminal will be restored a 
command will not be execute 



selection of a file on the user's 
the output file device. The 
and all arror messages will be 
ardless of the output file device 
equent ouput will be written into 
's A-disk. If N1 is 6, the 
s the output file device. This 
d on the MVS (batch) system. 



TIME 

This operation permits the time printout to be suppressed 
without loss of other output. A second TIME will restore 
the time printout unless the print output is suppressed with 
the NO command. 



SAVE or SAVE, Ml or SAVE,M1,N1 

fhis operation saves all arrays in storage at the time of 
issuance. Saved arrays will contain all modification^ made 
since their creation. Ml and N1 are optional and if not 
included, the arrays will be stored in FILE 02 on the user's 
A-disk. The saved files will be assumed to be for general 
matrix manipulation or a structural problem. Ml is the file 
name (up to six letters) under which the user wishes to 
store the arrays. The file type will be SAVE. If N1 is 1, 
a general matrix manipulation or structural problem is being 
saved. If N1 is 2, a heat transfer problem is being saved. 



RESUME or RESUME, M1 or RESUME, Ml, N1 

This operation reads a ?aved file into memory. Any arrays 
currently in storage will be destroyed. A file must have 
been previously created on the user's A-disk using the SAVE 
operation. Ml and N1 are optional and if not included, FILE 
02 will be read. It will be assumed to be a general matrix 
manipulation or structural problem. M1 is the file name 
assigned to the saved file on the user's A-disk. The file 
type must be SAVE. If N1 is 1 , a general matrix manipula- 
tion or structural problem is being resumed. If N1 is 2, a 
heat transfer problem is being resumed. 



LIST 

The LIST operation prints the directorv information for 
arrays in storage and the amount of storage used. 
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+ 

LOAD,m ,N1 ,N2,N3 

This operation will load an array of real numbers named Ml 
which has N1 rows and N2 columns. The terms of the array 
are punched in row-wise sequence on data cards following 
this operation. N3 is optional. If S3 is zero or blank, 
the cards are punched in (oP 10.0) format. If S3 is one, an 
additional card which contains the format of the data cards 
must preceed the data. Pop example, if the dat^ ps to be 4 
numbers per card in field widths of 15, the additional card 
would contain the following information: (4F15.0) . If S3 is 

nine, the data cards will Be read in free format. 



+ 



+ 



LOADI,M1,N1,N2,N3,N4 or LOA DI,M1 ,N1 ,N2 
LOADI,M1,N1,N2,N3 or LOA DI,M1 , SI ,N2 
L0ADI,M1 ,M2,N1,N2,N3, N4 



, ,N4 or 
or 



This operation will load an integer array named HI which has 
N1 rows and N2 columns.. The terms of the array are punched 
in row wise sequence on data cards following this operation. 
M2, S3 and S4 are optional. If S3 is zero or blank, the data 
must be punched in (1615) format. If S3 is one, an addi- 
tional card containing the format of the data cards must 
follow this operation and preceed the data. For example, if 
the data is to be 4 numbers per card in field widths of 10, 
the additional card would contain: (4110). If S3 is nine, 

the data will be read in free format. 

If the letters "ISCB” are placed in the position of M2, this 
operation has an increment generation capability. Data must 
be entered as follows: 

Item Contents 

1 Row number 

2 Value 1 

3 Value 2 

S2+1^’value S2 
S2+2 Generation code 

If the generation code is not zero, the next card must 
contain the following: 

Item Contents 

1 Row number increment 

2 Value 1 increment 

3 Value 2 increment 

0iic • 

S2+1 Value S2 increment 
S2+2 Last row to be generated 



+ 

ZERO,M1,S1,S2,S3,S4 



A real matrix named Ml is 
columns. The terms in this 
values: 



Therefore 

matrices. 



Ml (1,1) = S3 I = 

Ml (I, J) = S4 J = 

this operation can 



created 

matrix 



1 

1 



fir 



.SI 

.S2 

used 



with SI 
will have 



to form 



rows and N2 
the following 



null or unit 
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PRINT-MI or PEINT,M1,N1 or PRINT, Ml , N 1 , M2 or 

PRI NT,M1,N1,N2,N3 

This operation w4.ll print the array named M1 in a matrix 
format of up to exght columns per line. M1,N2,N3, and N4 
are optional. N1 is the number of comment cards (following 
the operation card) which will be read and printed. Nl 
defaults to zero. If N2 is included, the matrix will be 

E rinted in partitioned form with N2 columns per partition. 

ines will have N2*15 5 characters. N2 defaults to 8, 

printing 125 charcters per line. If M3 is included and is 
greater than zero, integer format (161 is used. The default 
value is zero and real format (PD15.7) is used. The user is 
cautioned not to overcome the capacity of the displaying 
device in use to avoid wrap around on the screen. 



♦ 

DUP,M1,M2 

Ttiis operation will form an array named M2 which is iden- 
tical to the array named Ml . 



ADD, HI, M2 

This operation will replace matrix Ml with the sum of the 
matrices Ml and M2. 



SOB, Ml, M2 

This operation will replace matrix Ml with matrix Ml less 
matrix M2. 



+ 

MOLT, HI , M2, M3 

This operation generates a new matrix M3 which is the 
product of matrices Ml and M2, or M3 = M1«^M2. 



+ 

TRAN, Ml, M2 

This operatign generates a new matrix M2 which is the trans- 
pose of matrix Ml. 



SCALE, Ml, M2 

This operation replaces each term in the matrix named Ml 
with the term multiplied by the term M2 (1,1) of the matrix 
named M2. 
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S0LVErM1,M2,N1,N2 or SOLVE, Ml ,M 2, , N2 or SOLVE, Ml , N1, N2 or 
SOLVE, Ml, M2, N1 or SOLVE, Ml, M2 

If Ml = 0, this operation solves the m^itrix equation AX=B. 
Ml IS the name of the A matrix and M2 is the name of the B 
matrix. Matrix A is triangularized and the results, X, are 
stored in M2. If N1=1, Matrix A is triangularized only. 

If N1=2, for a given B matrix and the A matrix previously 
triangularized, the B matrix is replaced by the results, X. 
If N1=3, Matrix A is replaced by its inverse FOR SYMMETRIC 
MATRICES ONLY. 

If N2=0 or blank, matrix A is symmetric. If M2 is nonzero 
the matrix A is not symmetric. 

For symmetric matrices, A is factored into the LDL form. 
The diagonal D matrix is stored on the diagonal of A. The 
parameter M2 permits the direct solution of non-symmetric 
systems of equations. If M2 is not equal to zero, an LO 
decomposition of matrix A will be performed. No direct 
replacement of Ml by its inverse is available for the non- 
symmetric case. Instead. use the ZERO operation to create 
an identity matrix M2 of the same order as Ml. The command 
SOLVE,M1 ,M2, . N2 will then replace the matrix M2 with the 
inverse of the matrix A. 



D0PSM,M1,M2,N1,N2,M3, N4 

This operation forms a new submafriy named M2 .with N3 rows 
and N4 columns from the terms within the matrix named M1. 
The first term of matrix M2, M2 (1,1). will be from row N1 
and column N2 of matrix Ml, or M1(N1 ,n2). 



ST0SM,M1,M2,N1,N2 

This operation s1:ores a submatrix named M2 within the matrix 
named HI. The first term of the submatrix M2 will be stored 
at row N1 and column N2 of matrix Ml. The terms within the 
area of Ml in which M2 is stored will be destroyed. 



OOP DG, Ml, M2 

This operation forms a new matrix named M2 from the diagonal 
terms of Ml. 



STODG,M1 ,M2 

This operation stores a row or column matrix named M2 at the 
diagonal locations of matrix Ml. 



+ 

MAX, Ml, M2 

This operation forms, a column matrix named M2 in which e^ich 
row contains the maximum absolute value of the corresponding 
row in matrix Ml. The maximum and its column number is 
printed for each row. 
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NORM, Ml, M2, N1 

If N1 = 0, a row matrix named M2 is formed in which each 

column contains the sum of the absolute values of the corre- 
sponding column of matrix M 1. If N1 is not egual to 0, a 
row matrix named M2 is formed in which each column contains 
the square root of the sum of the squares of the values of 
the corresponding columns of matrix HI. 



SQREL,M1 

This operation replaces each term in the matrix named Ml 
with the square root of the term. 



LOG , M 1 

This operation replaces each term in the matrix named Ml 
with the natural log of the term. 



PROD, Ml, M2 

This operation forms a 1 x 2 array named M2 which contains 
the product of all terms in the matrix named Ml. The 
product, X, is stored as two numbers of the form: 

X = psi«10S«^5E 

in which 112(1) = P and M2(2) = E, the exponent. 



DELETE, Ml 

This operation will cause the elimination from storage of 
the array named Ml. 



CMS ,N1 

This operation allows the user to issue CMS commands while 
under the control of CAL-NPS. In general a command that 
reloads the virtual core will not be allowed. Examples are 
FORTHX. SCRIPT, XEDIT, any language processor, SORT, and 
other large system modules such as COPiFILE and MOVFILE. N1 
is the number of words in the command (1 to 9 are allowed). 
Note that a parenthesis - or ”) ” is counted as a word. 
All words must be left justified. 

The user will be prompted for each word. If this operation 
IS used in an FT04F001 FILE, each word must be on a separate 
line. 
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• Stati c Analysis Ope rati o ns 

The purpose of this series of operations is to form 
the total stiffness and diagonal (lumped) mass matrices for 
systems of two- or three-dimensional elements. For three- 
dimensional analysis there are beam and truss elements 
available. For two-dimensional analysis, there is a frame 
element, a slope/deflection element for beams, and a 4 to 8 
node isoparametric finite element available. 

After the creation of an array containing the coor- 
dinates of the joints of the system, the specification of 
displacement boundary conditions, the tabulation of material 
and section properties, the mass and stiffness matrices are 
formed for each structural member and placed in sequence on 
low speed storage along with the global equation numbers 
which are associated with their stiffness terms. In addi- 
tion, the member force-displacement transformation matrices 
are formed and stored on a separate low speed storage file 
along with the appropriate displacement numbers. 

The NODES operation is used to specify or generate 
the geometry of the system. The operation BOUND specifies 
which joint displacements exist and assigns internal equa- 
tion numbers to these displacements. Therefore, each joint 
may have from zero to six displacement degrees of freedom. 
Tables of material and section properties for the various 
members are loaded and printed as standard arrays of 
information. 

A special operation, ADDSF, is used for the direct 
addition of element stiffnesses to form the total stiffness 
and diagonal mass matrix of the system. The ADDK operation 
may be used to add individual elements into the total system 
matrices. The LOADS operation specifies the concentrated 
joint loads for all load conditions. After the direct solu- 
tion for joint displacements due to static or dynamic loads. 



45 



the member forces can be evaluated using the MEMFRC opera 
tion. The DISPL operation is used to print the displace 
ments in joint number order. 
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+ + 

NODES, Ml, N1 or NODES, Ml, Nl, N2 

This operation generates an array (Nil, 3) .named Ml which 
contains the coordinates for all :)Oints in a structural 
system. N2 is optional. Da±a must be entered in free 

format as follows; 

Item Contents 

1 Node number 

2 X-coordinate 

3 Y-coordinate 

4 Z-coordinate 

If N2 = 1, there is a joint generation and coordinate system 
conversion capability. Data must be entered in free format 
as follows: 

Item Contents 

1 Node number 

2 x-coordinate 

3 Y-coordinate 

4 z-coordinate 

5 System type 

6 Generation code 

System type refers to the system used when inputing the 
data. All coordinates will be converted to the cartesian 
system for use by CAL-NPS. 

System type System 

“ Cartesian 

Cylindrical, Z axis longitudinal 
Cylindrical, Y axis longitudinal 
C^indrical, X axis longitudinal 
S pherical 

is the same as above with the following 



Cy lindr i cal Sp herical 



If the generation code is not zero, the next card is a 
generation vector for the self generation of nodes. It is 
formatted as follows; 

Item Contents 

1 Node number increment 

2 X increment 

3 Y increment 

4 Z increment 

5 Last node number to be generated. 

It is assumed that the increments pertain to the same system 
of reference as the preceeding card. 

This operation must be terminated by a line of alternating 
zeros and blanks. 



1 

2 

3 

4 

5 

The input data 
correspondence; 

Cartesian 

X 

Y 

Z 
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BOUND, Ml 

This operation specifies the displacements which are nonzero 
for the structural system of joints specified by the NODES 
operation. Where: 

Ml = Name of boundary condition code array to be 
gen erated. 

This operation is foliowetj by a series of cards containing 
the following information in free format: 

Item Contents 

1 Node number for the first node in a 
series of nodes with identical 
displacement specification. 

2 Nod$ number for the last node in the 
series. 

3 X-translat ion 

4 I-translat ion 

5 Z-translation 

6 X -rot at ion 

7 Y-rotation 

8 Z-rotation 

9 Node.ngmber increment used to generate 
conditions for addtional nodes. 

A translation or rotation equals: (a) zero foy zero or unde- 
fined displacements, or (b) one for nonzero displacements to 
be evaluated by other operations. 

If a node boundary condition is not specified, all displage- 
menrs at that node are assumed zero. Cards may be supplied 
in any order. If node boundary conditions are specified 
more than once, the last definition is used. This sequence 

of data must be terminated by a card of alternating zeros 
and blanks. 

Tbe selection by the user of which nodes h^ve nonzero 
displacments requires an understanding of the direct stiff- 
ness procedure. Displacement degrees of freedom which have 
no stiffness associated with the displacement must be 
considered to be undefined since it is not possible to 
develop an equilibrium equation for that direction. The 
total number of nonzero displacements specified will be the 
size of the total stiffness matrix to be defined by the 
ADDSF operation. 
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BEAM, Ml , M2, M3, H4 

This operation calculates the element stiffness, mass and 
force-aisplacement transformation matrices for 3-D beam 
members. These arrays are stored in sequence on low speed 
storage to be used by other operations where: 

Ml is the name of the beam element group 

M2 IS the name of the coordinate array 

M3 is the name of the boundary condition array 

M4 is the name of the array which contains beam 

properties and has been loaded by the standard 
matrix LOAD operation 

One card for each. beam in this group of beam elements must 
follow this operation. The Beam cards are punched in free 
format, where: 



Item Contents 

1 Beam identification number 

2 Node number I 

3 Node number J 

4 Node number K 

5 Beam property number NP 



This sequence of cards, must be terminated with a card of 
alternating zeros and blanks. 



The material and geometric properties for each element are 
given in the M4 array in the following order: 



H4 

M4 

M4 

M4 

M4 

M4 

M4 



(NP, 1) 
/NP,2 
’np' 3 
NP,4) 
NP,5 
NP,6 
NP,7! 



Axial area of member, A 

Torsional moment of intertia, J 

Moment of inertia about axis 2, I 

Moment of inertia about axis 3, I 

Modulus of elasticity, S 

Shear modulus, G 

Mass per unit length of beam 



whepe NP is the specific material property number spegifi^d 
in Item 5 of the beam card. The local sign convention is 
given in the following figure. 
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2 






Ya 




X 





LOCAL BEAM REFERENCE SYSTEM: 
AXES 1 AW 2 ARE IN THE 
PLANE DEFINED BY PLANE 
I J AND K NODES. AXIS 1 
IS DEFINED BY LINE I-J. 
AXIS 3 IS PERPEWICULAR 
TO THE 1 , J , R plane. 



LOCAL SIGN CONVENTION 




DEFINITION OF POSITIVE BEAM FORCES 
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♦ 

TBOSS^fll ,a2rH3r U4 



This operation forms the. element. stiffness, mass and force- 
displacement transformation matrices for 3-D truss members. 
The arrays are stored on low speed storage in sequence and 
will be used by other structural operations. 

HI is the name of this group of truss members 
H2 IS the name of the coordinate array 
H3 is the name of the boundary . condition array 
HU is an HP by 3 array of section properties in 
which HP is the number of different section 
properties and 



H4(NP,1) a The cross-sectional area, A 
H4(NP,2i a The modulus of elasticity, E 
H4 (HP, 3) a the mass per unit length of the member. 
This matrix can be loaded by the matrix LOAD operation. 



This 

free 



operatioi} is followed. by <?ne card, per truss member in 
format with the following information: 

Item Contents 

1 Truss member identification number 

2 Joint number I 

3 Joint number J 

4 Section property number, HP 



This operation must be terminated by a card of alternating 
zeros and blanks. 



PRAHE,H1,M2 

This opera'tion forms the 6 x 6 stiffness matrix for the 
two-dimensional frame member shown below. 




FRAME MEMBER 
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The properties of the neaber are defined on one cards xaae- 
diately following the FRAME operation card. This second 
C 2 ird is punched in free format and contains the following 
inf or nation: 

Ite* Coisitents 

1 Azxal area, A 

2 Modulus of elasticity, E 

3 Moaent of inertia, I 

4 Xi 

5 Yi 

6 ij 

7 Ti 




M2 is a 3 X 6 force-displac eaent transformation matrix which 
is based on the positive definition of the element forces 
shown below. 
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These forces can be calcu lated froe 
equation with the flEHFBC operation. 

ru*7 



M 

M 

P 



= [M2] 



Uyi 

U»4 

U.| 

Uy| 

U«i 



the following matrix 



SLOPE, Ml 

This operation forms a u x 4 stiffness matrix, B1 fpr a beam 
or coluon member from the classical slope'-deflection equa-- 
tions. The properties of the member are defined on one card 
immediately following the operation. This second card is 

f unched in free format and contains the following informa-* 
ion: 

Item Contents 

1 Homent of inertia, I 

2 Hodulus of elasticity, E 

3 Length of member, L 

The sign convention is defined as follows: 





DISPLACEMENTS 



The member forces are defined in terms of joint displace- 
ments by the following slope deflection equations. 



II 




M,= il' 




II 

1 

II 





L 
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PLANE^HI M4,N1, N2 

This operation calculates the element stiffness, mass and 
stress-displacement transformation matrices for 4 to 8 node 
isoparametric elements (Y-Z plane only). -These arrays are 
stored on low speed storage to be used later by other opera- 
tions (i.e. , ADDSF and FORCE). The arguments are defined as 

Ml is the user defined n^me of the element group 

M2 is the name of the joint coordinate array 

M3 is the name of the boundary condition array 

M4 is the name of the array which contains the 

material properties of the elements (one row per 
different material) where 



M4 

M4 

M4 

M4 



(NP, 1 
'nP,2 
NP,3 
NP,4; 



Modulus of elasticity, E 
Poissons ratio, v 
Thickness of element 
Mass density of the element. 



NP is the material identification number. 



N1 5 ind N? are the nupber of integration points in the r and 
s directions respectively. 

One card for each 3 to 8 node element in the gcoijp must 
follow the operation card. The cards are punched ifi free 
format and contain the following information: 

Item Contents, . . 

1 Element identification number 

2 Node number N1 

3 Node number H 2 

4 Node number N3 

5 Node number N4 

6 Node number N5 

7 Node number N6 

8 Node number N7 

9 Node number N8 

10 Material identification number, NP 

11 Natural Coordinate of stress output r1 

12 Natural Coordinate of stress output si 

13 Natural Coordinate of stress output r2 

14 Natural Coordinate of stress output s2 

15 Natural Coordinate of stress output r3 

16 Natural Coordinate of stress output s3 

N4 through N8 are optional, but zeros must be inserted for 
them, if unused. The midside nodes, if present, must be 
within the center half of the side. The local numbering 
system for the element is shown in the following figure. 
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ISOPARAMETRIC ELEMENT 

Stresses will be printed by the FORCE operation at rhe three 

J oints defined in iteas 1 1 through 16. The forces are 
efined as follows: 
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+ 

LOADS, Ml, M2, N1 

This operatioij .forms a load matrix named Ml of N1 columns 
(N1 load conditions) where M2 is the name of the boundary 
condition array generated by the operation BOUND. This 
operation is followed by a series of cards - one for each 
loaded joint for each load condition. These cards are 
punched in free format as follows: 



Item Column 

1 Joint number 

2 Load condition number 

3 Load in X-direction 

4 Load in Y-direction 

5 Load in Z-direction 

6 Momemt about X-axis 

7 Momemt about Y-axis 

8 Momemt about Z-axis 



This series of cards must be terminated by row of alter- 
nating zeroes and blanks. 



♦ + + 

ADDSF,M1 or ADDSF,M1,M2 

This operation forms the total stiffness matrix named Ml and 
the lumped mass matrix named M2 for the structural system 
from the element stiffness and mass matrices which are 
stored on low speed storage. These matrices can be printed 
with the PRINT operation. If M2 is not specified, the row 
mass matrix M2 will net be formed. 



ADDK,H1 ,M2,M3,N1 

This operation adds the element stiffness matrix named M2 to 
the total stiffness matrix named Ml, where Ml was previously 
defined and initially set to zero, M3 is the name of the 
integer array in which the column number N1 contains the row 
or column numbers in the total stiffness matrix where the 
element stiffness terms are to be added. 



DISPL,M1,M2 

This operation prints the displacement vector named Ml in 
joint sequence order, where M2 is the name of the boundary 
condition array. 
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MEHFEC^MI, M2, M3,M4, N1 

This operation multiplies the element stiffness matrix named 
Ml by the joint displacement matrix named M2. M3 is the 
name of the integer array in which the column number N1 
contains the row numbers in the displacement matrix, M2, 
which are to be multiplied by the element stiffness (or 
force-displacement) matrix. Ml. The results of this multi- 
plication are stored in the array named M4. 



FORCE, Ml , M2, M3 or FORCE, Ml, M2 

This operation calculates the member forces for a group of 
elements in which 

M1 is the name of the element group 

M2 is the name of the displacement matrix 

m 3 is the name of the matrix in which the forces are 

stored in the order calculated. 

If M3 is not specified, tfeie element forces will be printed 
only and will not be retained in storage. For the TROSS 
element only the member axial force, F, will be calculated 
for each member. For the BEAM element, eight forces will be 
printed with reference to the positive definition shown in 
the BEAM operation. 
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3. Dynamic An alysis Op eration s 

The following operations were designed to eyaluate 
the dynamic response of structures subjected to arbitrary 
time- dependent loads. If these operations are used in 
connection with the standard matrix operations and the 
structural analysis operations, a dynamic analysis is a 
relatiyely simple procedure. The user has the option of 
using the mode superposition method or a direct step-by-step 
integration of the dynamic equations of motion. The user 
may examine the spectra of both input loading and calculated 
displacements. In addition, the contributions of the indi- 
yidual modes may be evaluated and compared. 

The most common and convenient form for time- 
dependent data to be specified is as straight line segments 
between given time points. Therefore, an operation which 
generates values at equal intervals is necessary. Another 
common characteristic of time-varying loads on structures is 
that it is normally possible to represent the loads at all 
points on the structure by the product of two matrices, a 
column matrix indicating the spatial distribution of loads 
times a row matrix which indicates the values as a function 
of various times. If a more complicated loading is 
required, it is possible to perform more analyses, each 
within the restrictions of the program, then add the results 
of each analysis. 

The following operations have been added for the 
major purpose of performing dynamic analysis. 
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♦ 

FUN6,ai,l!2,I13,H1,ll2 

This operation generates a natrix named M2 which contains 
values at equal intervals of the function specified in the 
array named Ml. The array ai must be a 2 by K array of the 
fora: 




which numerically represents a function of the form shown 
belov: 



F(T) 




The tine interval t is specified in the 1 by 1 matrix named 
M3* HI specifies the total number of values to be gener- 
ated, and is the number of columns in M2. If N2 = 0, the 
array M2 will be a 1 x HI row matrix in which the first 
value will be f. If N2 is not equal to 0, the array M2 will 
be a 2 X HI natrix of the following form; 



[M2] 



t, t, ♦ At ti ♦ 2At 
f, f (t, ♦ At) f (t, + 2At) 
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STEP^MI , M2, 113,114, M5,M6,M7, M 8, N1 , N2 

This operation calculates the dynamic response of a struc- 
tural system using direct step-by-step integration of the 
following linear matrix equation of motion: 

[M]{a} + [C]{0} ♦ CK]{U) = R(t) = PF(t) 



Hhere : 



Ml is the name 
M2 is the name 
M3 is the name 
H4 is the name 
D in which: 




of 

of 

of 

of 



the 
t he 
the 
the 



V ect or 

V ect or 
vector 



N 

N 

N 

N 

Of 

of 

Of 



X N stiffness matrix K 
X N mass matrix M 
X N damping matrix C 
X 3 initial condition matrix 



disulacements 
velocities U 
accelerations 



0 



M5 i$ the name of the by N2 matrix of calculated 
displacements in which column i represents the 
displacements at time i=^N1* t 
M6 is the name of the N x 1 load distribution matrix 
P 

M7 is the name of the 1 x k row matrix representing 
the load multipliers at equal time increments 
P, where k = N2/N1 

M8 is the name of the 1x1 matrix containing t 
N1 is the output interval for the displacements 
N2 is the total number of displacement vectors to 
be calculated. 



The total tioe for which results will be calculated by this 
operation is N1’i‘N2‘^ t. This operation must be followea with 
one data card in free format containing the following infor- 
mation: 



Item Contents 

1 DELTA 

2 ALPHA 

3 THETA 



Different values gf delta, alpha and theta will a^low the 
user to select different methods of step-by-step integra- 
tion. The following table lists some possibilities: 



DELTA ALPHA THETA 



New marks Average Acce 
Linear Acceleration 
Wilson's Theta Method 
Wilson's Theta Method 



lerati on 

flow damping) 
(high damping) 



1/2 

1/2 

1/2 

1/2 



1/4 

1/6 

1/6 

1/6 



1.0 

1.0 

1.42 

2.0 
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EIGEN,Ml,M2,M3rN1 



This operation solves the following eigenvalue problem: 

K $ = M 4) A 



In which the N x, N, symftri 
IS named Ml. The matrix H 
positive terms designated 
row or column matrix conta 
M. The eigenvalues, , are 
values are ordered in nume 
eigenvectors, , are store 
the matrix M2. The numb 
number of significant figur 
zero or blanJc, 4 figure acc 
accuracy possible is 16 fi 
figure accuracy is not reco 



c, positive-semidefinit § matrix K 
IS a diagonal matrix or nonzero, 
by M3. The matrix M3 must be a 
ining only the diagonal terms of 
stored in matrix m 3. The eigen- 
rically increasing order and the 
d in the corresponding columns of 
er N1 specifies the approximate 
es of the eigenvalues. If N1 is 
uracy will be used. The maximum 
gures. The use of more than 12 
amended. 



The program reduces the problem 
by the following transformation 

where 



to standard eigenvalue form 



m = m K m 
I = m ^ M m 

in which 

m, = 1/ /Mil 

The calculated mode shapes, , are normalized as follows: 

The program uses the standard Jacobi diagonalization method 
to solve for all eigenvalues and eigenvectors. 
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+ 

DYNAM,tl1,il2,I13, M4,M5, M6, N1 



This operation ^evaluates 
second order differential e 
superposition method for t 
tural system. 



the .following, set of uncoupled 
quations associated with the mode 
he dynamic analysis of a struc- 



X, + 2A,tD,x, + wfx, = F?(t) i = 1 to N nodes' 



Ml is the name of a row or column matrix which contains the 
N terms (frequencies in rad/sec) . M2 is the name of a row 
or column matrix which contains the N terms (ratio of 
modal damping to critical damping). 

The generalized time-varying forces P (t) are not specified 
directly but are evaluated from more fundamental informa- 
tion. The forces for all modes are evaluated at specific 
times by the program from the following matrix equation: 

p = p^f = M3*'H4 

In which p is a specified {1x1 vector named M3, and f is a 
1 X N1 row matrix which will be generated from the 2 by k 
array named H4. The array M4 is the same form as the input 
array described under the operation FONG. It is not neces- 
sary to use FONG before the DYNAH operation. 

M5 is the name of the N x N1 array which contains the gener- 
alized displacement X (t) . 

M6 is the name of the. 1 x 1 array which contains the time 
increment associated with the generalized displacements. 

N1 is the number of displacements to be generated. 

The method of integration used is exact for straight line 
segments. 



PLOT, Ml, N1 

This operation will prepare a pointer plot of selective rows 
of the matrix named Ml. N1 is the number of rows of Ml 
which will be plotted by this operation. This operation is 
followed by N1 cards in f (1a1,i4) format with the following 
information: 

Columns Contain 

1 Plot symbol - any keypunch symbol 
2 - 5 Bow number to be plotted 

The program automatically searches the information to be 
plotted for the maximum and minimum values. The difference 
in these numbers divided by 120 spaces is selected as the 
plot scale. 
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^ Transfer Oper a tion s 



The purpose of this series of operations is to form 
the total conductivity and heat capacity matrices for 
systems of two- or three-dimensional elements, to form the 
flux vector and solve the defined set of equations. For two 
dimensional elements there are 4 to 9 node isoparametric 
elements. For three dimensional elements there are 8 to 21 
node isoparamtric elements. 

After the creation of an array containing the coor- 
dinates of the system nodes, the specification of element 
connectivity, the specification of material properties and 
the specification of constant temperature nodes, the conduc- 
tivity and capacity matrices and the flux vector are formed. 
The equations are solved by the appropriate equation solver 
and the temperatures are printed in node order. 

The HTXFR operation initializes the problem. The 
operation COORD is used to specify or generate nodal coordi- 
nates. Element connectivity is specified by the ELCON oper- 
ation. Material properties and element boundary conditions 
are input via the PROP operation. The CTEMP operation 
establishes designated nodes as having constant tempera- 
tures. The equation profile for the problem is generated by 
the PROF operation. 

To form the conductivity matrix, the operation SYHC 
is used. The heat capacity matrix can be approximated with 
either a consistant (CCAP) or lumped (LCAP) matrix formula- 
tion. The flux vector is formed with the FORM operation. 

Time independent systems of equations are solved 
with the CALC operation. Systems of equations involving the 
first derivative of temperature with respect to time are 
solved with the ODE series of operations. Nodal temepra- 
tures are printed in node number order with the PTEHP opera- 
tion. 
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HTX FR 



? his operation initializes the h$at , transfer problem., , Jbe 
ollowing information, entered in free format must follow 
this operation: 

Item Contents , 

1 Number of nodes (NOMNDP) 

2 Number of elements (NUML) 

3 Number of material sets (NOMMAT) 

4 Spatial dimension (2 or 3) (NDMi 
5 Number of unknowns per node (NDF) 

6 Maximum number of nodes per element (NEN) 

The number of unknowns per node will always be 1 for heat 
transfer problems. 



COORD 

This operation creates an array which contains the coordi- 
nates for all nodes in a heat transfer system. Data is 
entered in free format as follows: 



Item Contents 

1 Node number 

2 X-coordinate 

3 Y- coordinate 

4 Z-coordinate (if 3-D, else omit) 

5 System type 

6 Generation code 

System type refers to the coordinate system used when input- 
ting data. All coordinates are converted to the cartesian 
coordinate system for use by CAL-NPS. 



System Type System 

0,1 Cartesian 

2 Cylindrical, Z-axis longitudinal 

3 Cylindrical, Y-axis longitudinal 

4 Cylindrical, X-axis longitudinal 

5 Spherical 



The input data is the same as shown above with the following 
correspondence: 



Cartesian 

X 

Y 

z 



Cy lindrica 1 
r 
0 
Z 



Soherical 

P 

0 

<t> 



If the generation code is not zero (0) , the next card is a 
generation vector for automatic node generation. It is in 
free format as follows: 



Item Contents 

1 Node number increment 

2 X increment 

3 Y increment 

4 Z increment (if 3-D, else omit) 

5 Last node number to be generated 

This operation must be terminated by a card with alternating 
zeroes and blanks. 
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ELCON 



This operation creates an. array which <;ontains the element 
connectivity data for the elements in a heat transfer 
system. Data is entered in free format as follows; 



Item 


Contents 


1 


Element number 


2 


Node 1 number 


3 


Node 2 number 


etc 


etc 


N-M 


Node N number 


N*2 


Material set number 


N+3 


Generation code 



If the generation code is not zero. (0) , the. next card is a 
generation vector of the automatic generation of element 
connectivity. It, is entered in free format as follows: 

Item Contents 

1 Element number increment 

2 Node 1 increment 

3 Node 2 increment 

Qt C 

N-*-l Node N increment 

N+2 Haterial set increment (usually 0) 

N+3 Last element number to be generated 

This operation must be terminated by a row of alternating 
zeroes and blanks. The node numbering conventions for 
element connectivity are shown below. 



4 7 5 

T 



8 9* t>S 6 



1— 5 2 




2-D Element 



3-D Element 
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PROP 



This operation inputs the material property data for. a heat 
transfer system. The following information, entered in free 
format, must follow this operation: 

Item Contents 

1 Material set number 

2 Element type number (2 for 2-D, 3 for 3-D) 

additional information must be provided, depending on the 
type of element being used. Note that material property 
data must be in consistent unit s as shown below: 



k 

c 

rho 

q 1 1 I 



English Units 

BTU/lbm-op 

lbm/ft3 

BT0/hr-ft3 



SI Uni ts 
— w7Tr=^c 
kJ/kg-oc 

kg/m3 

w/m3 



2-D Eleme n ts 

This information must follow the PROP operation card and its 
two required entries. An entry must be made for each item. 
If the item is temperature dependent, the entry will be 
ignored. 

Item Contents . 

1 Conductivity in the X-direction 

2 Conductivity in the Y-direction 

3 Specific Heat 

4 Density 

5 Heat generation per unit volume 

6 Number of integration points per direction 
(1 to 6, default is 4) 

7 Geometry type (see below) 

8 Total number of lines with specified boundary 
conditions in elements with the same material 
set numer (NLBC) (see below) 

9 Temperature dependence code (see below) 

Geometry type is 1 for plane geometry and 2 for axisymmetry. 

The temperature dependence gode is 0 if all material proper- 
ties are constant and 1 if any property is temperature 
dependent. If the code is 1, the following information (in 
free format) is required. 

Item Content? 

1 Conductivity in the X-direction code 

2 Conductivity in the Y-direction code 

3 Heat capacity (specific heat* density) code 

4 Heat generated per unit volume code 

where 0 means a constant property and 1 means a temperature 
dependent property. Temperature dependent properties are 
entered in the form of a table. The tables are consecu- 
tively input for conductivity in the X-direction, conduc- 
tivity in the Y-direction, heat capacity and heat generated 
per unit volume. Omit the tables for which the temperature 
dependence code is zero. Tables are input in free format as 
shown on the following page. 
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I tea Contents 

1 Humber of data pairs to be entered (This 

should be a single card) 

2 Temperature 1 

3 Property 1 (These two entries should be on 

one card) 

4 Temperature 2 

5 Property 2 (These two entries should be on 

one card) 

0*tc 0iic 

2*H Temperature H 

2<'H+l Property H (These two entries should be on 

one card) 



If any lines have, a specified boundary condition (MLBC;>0),,a 
card must be submitted for each line. If the same line is 
subjected to more than one boundary condition, a card must 
be used for each one of these conditions- The total number 
of cards must equal HLBC. The information is entered in 
free format as follows: 

Item Contents 

1 Element number 

2 Boundary condition code (see below) 

3 Line code (see below) 

4 Property value (see below) 

5 &mb;.ent temperature 
The boundary condition codes are: 

1 Flux 

2 Convection (constant coefficient) 

3 Radiation 

4 Convection (temperature dependent property) 

The line codes are: 

1 S * -M line 

-1 S s -1 line 

2 T = Vi line 

-2 T * -1 line 



4 7 3 

T 

4 



8 9 oS Q 



5 2 



The property values are: 

Flux - Flux per unit area 

Convection - constant heat transfer coefficient 
(ignored if temperature dependent) 

Radiation - Product of emissivity by Stephan- 
Boltzman constant 

The ambient temperature is ignored for the flux boundary 
condition. 

If the boundary condition code is 4 (temperature dependent 
heat transfer coefficient) - a table musr follow for the 
temperature dependence- It is input in free format as is 
shown on the following page. 
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Item Contents 

1 Number of data pairs to be entered (This 
should be a single card) 

2 Temperature 1 . . 

3 Heat transfer coeficient 1 (These two entries 
should be on one card) 

4 Temperature 2 

5 Heat transfer coeficient 2 (These two entries 
should be on one card) 

©t c © t c 

2^N Temperature H . 

2*N+1 Heat transfer coeficient N (These two entries 
should be on one card) 



3-D Eleme nts 

This infgrmation ipust follow the PROP operation card an4 its 
two required entries. An entry must be made for each item. 
If the item is temperature dependent, the entry will be 
ignored. 

Item Content^ . . , 

1 Conductivity in the X-direction 

2 Conductivity in the Y-direction 

3 Conductivity in the Z-direction 

4 Specific Heat 

5 Density 

6 Heat generation per unit volume 

7 Humber of integration points per direction 
(1 to 6, default is 4) 

8 Geometry type (see below) 

9 Total number of surfaces with specified 
boundary conditions in elements with the same 
material set numer (NSBC) (see below) 

10 Temperature dependence code (see below) 

Geometry type is 1 for plane geometry and 2 for axisymmetry. 

The temperature dependence code is 0 if all material proper- 
ties are constant and 1 if any property is temperature 
dependent. If the code is 1, the following information (in 
free format) is required. 



Item Contents 

1 Conductivity in the X-direction code 

2 Conductivity in the Y-direction code 

3 Conductivity in the Z-direction code 

4 Heat capacity (specific heat^* density) code 

5 Heat generated per unit volume code 



where 0 means a constant property and 1 means a temperature 
dependent property. Temperature dependent properties are 
entered in the form of a table. The tables are consecu- 
tively input for conductivity in the X-direction, conduc- 
tivity in the Y-direction, conductivity in the Z-direction, 
heat capacity and heat generated per unit volume. Omit the 
tables for which the temperature dependence code is zero. 
Tables are input in free format as is shown on the following 
page. 
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Item Contents 

1 Sanber of data paxrs to be entered (This 

should be a single card) 

2 Teaperature 1 

3 Property 1 (These two entries should be on 
one card) 

4 Temperature 2 . , , . 

5 Property 2 (These two entries should be on 
one card) 

0*t C 

2*H Teaperature M 

2*Nyl Property H (These two entries should be on 
one card) , . . . 

If any surfaces have a specified boundary condition 
(HSBOO) . a card must be su bait ted for each surface. If the 
same surface is subjected to more than one boundary condi- 
tion. a card must be used for each one of these conditions. 
The rotal nuaber of cards lUSt equal NSBC. The information 
is entered in free format as follows: 

Item Contents 

1 Eleaent number 

2 Boundary condition code (see below) 

3 Surface code (see below) 

4 Property value (see below) 

5 &ab;.ent teaperature 
The boundary condition codes are: 

1 Flux 

2 Convection (constant coefficient) 

3 Badiation 

4 Convection (teaperature dependent property) 

The line codes are: . 

1 a s +1 line 

-1 a = -1 line 

2 S « *1 line 

-2 S = -1 line 

3 T * +1 line 

-3 T * -1 line 




The property values are; 

Flux - Flux per unit area 

Convection - constant heat transfer coefficient 
(ignored if temperature dependent) 

Badiation - Product of emissivity by Stephan- 
Boltzaan constant 

The ambient temperature is ignored for the flux boundary 
condition. 

If the boundary condition code is 4 (temperature dependent 
heat transfer coefficient) . a table must follow for the 
temperature dependence. It is input in free format as is 
shown on the following page . 
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Item Contents 

1 Number of data pairs to be entered (This 
should be a single card) 

2 Temperature 1 

3 Heat transfer coeficient 1 (These two entries 

should be on one card) 

4 Temperature 2 

5 Heat transfer coeficient 2 (These two entries 

should be on one card) 

C 0tC 

2^N Temperature N 

2*^N + 1 Heat transfer coeficient N (These two entries 

should be on one card) 



CTEHP 

This operation inputs constant temperature boundary 
restraint data. If you have nodes with constant tempera- 
ture, you must enter the temperatures for those nodes. For 
nodes with no temperature restrictions, no entries should be 
made. Automatic generation capability is built into the 
operation. Data is entered in free format as follows: 

Item Contents 

1 Initial node 

2 Last node 

3 Node increment 

4 Temperature 

These entries may be repeated until al+ the constant temper- 
ature nodes are entei^ed. This operation must be terminated 
by a row of alternating zeroes and blanks. 



PROF 

This operation establishes the profile of the equations for 
solution of the heat transfer problem. After the issuance 
of this command, the problem is set and you may not change 
the node numbers with restrained boundary conditions 
(constant temperatures). The values of the restrained temp- 
eratures may be changed. 



SYMC 

This operation forms the symmetric conductance matrix for 
heat transfer problems. 



OSYMC 

This operation forms an unsymmetric conductance matrix. 



ICAP 

This operation forms a lumped capacitance approximation 
matrix for heat transfer problems. 
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ccap 



This operation forms a consistent capacitance approximation 
matri'X for heat transfer problems. 



FORM 

This operation forms the flux vector for heat transfer prob 
lems. 



CALC 

This operation solves time independent heat transfer prob 
lems for temperature and updates the temperature matrix. 



0DE,M1 

This operation. solves the. first order ordinary differential 
equations arising from time dependent heat transfer prob- 
lems. If Ml is ”INIT”,the initial conditions are estab- 
lished. This operation must be followed by the following 
data in free format: 

Item Contents 

1 Integration paramter theta for two point 
scheme (default = 2/3) 

2 Integration parameter gamma for three point 
scheme (default = 1.5) 

3 Integration parameter beta for three point 
scheme (default = 0.8) 

4 Maximum temperature difference for time step 
adjustment 

5 Minimum temperature difference for time step 
adjustment 

Default values are obtained by entering zeros for items 1, 2 
and 3 . 



If Ml is "LINE”, the two point integratio^i scheme is used. 

If Ml IS '•QUAD”, the three point integration scheme is used. 



Some suggested values for theta for the two-point scheme as 
given in Reference 4 are: 



Crank- Nico Ison 
Zienkiewicz 
Bettencourt 
Liniger 



THETA 

1/2 

2/3 

3/4 

.878 



Some suggested values for beta and gamma for the three point 
scheme as given in Reference 4 are: 



Lees 

Hogge 

Hood 

Zienkiewicz 

Bettencourt 



Beta 

V3 

3/4 

.646 

4/5 

9/10 



Gamma 

1/2 

1.0 

1.2184 

3/2 

3/2 
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PTE HP 

This operation prints the nodal temperatures of a heat 
transfer system in node number order. 



T0L,H1 

This operation $ets the solution convergence tolerance to 
the value found in the 1x1 matrix named Hi. 



CONV 

This operation performs a temperature convergence test on a 
heat transfer system. If this operation is used inside a 
loop (LOOP operation) and the test shows convergence, 
looping will be terminated. 



DTI H, HI 

This operation sets the time increment for integration in a 
heat transfer system to the value found in the 1x1 matrix 
HI. 



ADTIH 

This oper?ition advances the time in ^ h$at transfer problem 
by one time step. The time step is input with the DTIH 
operation. 



EIGV 

This operation computes the dominant eigenvalue and eigen- 
vector of the current heat transfer conductance matrix. 



PBOHPT 

This operation permits the prompts for user input to be 
suppressed without loss of other output. A second PROHPT 
will restore the prompts unless the print output is 
suppressed by the NO command. 
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5 • Graphics Operations 



CAL-NPS has a limited .graphics capability for users 
with PLOT- 10 compatible terminals or IBH 3277 Dual Screen 
terminals. Two- or three-dimensional meshes may be viewed 
in the X-T, T-Z or X-Z planes. 

The graphics operations are initialized with the 
GRAPH operation by which the user specifies which graphics 
capable terminal is being used. The PLHX operation plots 
heat transfer meshes and the PLST operation plots structural 
meshes. The user nay title the plots using the TITLE opera- 
tion. 



GRAPH 

This operation initializes the graphics . package. Presently 
there are two terminals users may utilize for the graphics 
operations. A yes must be input immediately following this 
card to signify that one of the two types of terminals is 
being used. Then the terminal code must be entered in free 
format; 

Terminal Code Terminal 

1 PLOT-1 0 Compatible Terminal 

2 IBH 327 7 Dual Screen 



PLHX,N1 

This operation plots ,2-D and 3-D heat transfer analysis 
meshes. If the mesh is 3-D, the user may specify the plane 
of view; 

m = 1 ; X-T plane 

N1 = 2 ; Y-Z plane 

N1 = 3 : X-Z Plane 

The default is the X-T plane, hence SI is not required for 
2-D meshes. 
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PLST,i!1,a2,»1 



Till? operation plots, 2-D an 
til IS the node coordinate a 
tion. 82 is the eleient 
the LOiDI operation. Th« 
eleaents and the colaan dx 
nodes per element plus on 
information: 

l!2(K, 1) s Number of no 
82 Ik, 2) = Node 1 
82 (K, 3) = Node 2 

82{K^N+1) *^Node N 
N1 specifies the plane of 7 
N1 s 1 : X-i plane 

N1 = 2 : I-Z plane 

N1 s 3 : X-Z plane 

If the structure contains 
(PLANE operation) , the con 
tion shown below: 

Z 



d 3-D structural analysis meshes, 
atrix cheated by the NODES opera- 
connectiyity matrix- loaded with 
row dimension is the number of 
mension is the maximum number of 
e. 82 contains the following 

des in element number K 



lew: 



two dimensional membrane elements 
nectiyity must follow the conven- 




Y 



TITLE,N1 

This operation allows the user to input N1 (up to three) 
fifteen character lines to label plots generated with either 
the PLHX or PLST operations. The label will appear in the 
upper right hand area of the screen outside the plotting 
area. This operation must immediately preceed the plotting 
operation- The title may be changed by reissuing the TITLE 
command. The default value of N1 is one. 
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6 . Looping Operations 



CAL-NPS has a five level looping ability. The first 
operation is LOOP and the last operation is NEXT. 

Operations within CAL-NPS are normally executed as they are 
encountered. If the operation requires data, the data cards 
follow the operation card. In the case of looping, however, 
all operation cards are stored within the computer before 
they are executed. If operations within the loop require 
data, the data cards must be supplied in the order required 
after the last NEXT operation. If an error is encountered 
while executing in a loop, the entire matrix of loop 
commands is deleted and the user is given the opportunity to 
try again. Matrices that have been modified by operations 
successfully completed while in the loop remain modified. 
After all loops are executed the computer storage required 
for these operations is automatically released by the 
program. 
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Loop,m 

N1 is the number of times the loop is to be executed. 
Associated with each LOOP operation there must be a corre- 
sponding NEXT operation which signifies the end of the loop 
and the return of the control to the beginning of the loop. 
The following is a possible series of looping operations. 



LOOP, 5^ 



LOOP, 2- 



NEXT 



Second level loop 
executed twice 
(total of 10 times) 



LOOP, 4 

NEXT 
NEXT ■ 






Second level loop 
executed 4 times 
(total 20 times) 



First level loop 
executed 5 times 



NEXT, Ml or NEXT 

The operatipn NEXT signifies the end of a loop. . It is 
apparent which LOOP and NEXT cards are associated if there 
are an egual number of each. The operation NEXT, Ml will 
cause the loop to terminate if the first term in the matrix 
Ml is negative. 



SKIP, Ml, N1 

This operation will cause the. skip of the. next N 1 . operations 
if the first term in the matrix named Ml is negative. 



7. Dser Defined Operat ions 



DSERA and OSERB 

These names are reserved for operations to be defined and 
programmed by the user. In order to program these opera- 
tions it is necessary to understand the internal organiza- 
tion of CAL-NPS. chapter III of Reference 2 contains the 
det ails. 
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E. LARGE PBOBLEHS 



CAL-NPS is designed as an educational tool. It does not 
take advantage of banding and symmetry in matrix storage, 
except in the heat transfer operations. Larger problems can 
be solved by increasing the dimension of the L array, but a 
general purpose program that makes maximum advantage of 
out-of-core storage and takes advantage of banding and 
symmetry for in-core matrix storage is probably a better 
choice. With the above disclaimer, to increase problem size 
capability, increase the dimension of the L array and change 
the value of MAX to»the new dimension size in the following: 

C BAIN PROGRAB 

C SET PROGRAB CAPACITY 

COBBON BTOT,NDP, L (100000) 

BTOT = 100 000 
NDP = 2 
CALL SETIBE 
CALL CALI 
STOP 
END 

With the dimension of the L array as above, the program 
currently executes in 1024 K bytes for CP/CBS. The region 
necessary for execution will increase about eight times the 
increase in the L array. 
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APPENDIX B 



SAMPLE DATA FILE 

This is a sample data file (FILE FT04F001) for the 
hollow cylinder with circumferential heating strips problem. 



HTXFR 
99 80 1 
2 1 4 



COO 


RD 






1 , 


1666667 5.819103 


2 


1 


1 1 


-.0013019 0 23 






2 . 


1666667 4.252421 


2 


1 


11 


-.0013019 0 24 






3 . 


1666667 3.469080 


2 


1 


11 


-.0013019 0 25 






4 . 


1666667 2.909551 


2 


1 


11 


-.0013019 0 26 






5 , 


1666667 2,461928 


2 


1 


11 


-.0013019 0 27 






6 . 


1666667 2.238116 


2 


1 


11 


-,0013019 0 28 






7 . 


1666667 2.014305 


2 


1 


11 


-.0013019 0 29 






8 . 


1666667 1.790493 


2 


1 


11 


-.0013019 0 30 






9 . 


1666667 1.342870 


2 


1 


1 1 


-.0013019 0 31 






10 


.1666667 0.78334 


1 2 


11 


-.0013019 0 32 






11 


.1666667 0 2 1 






11 


-.0013019 0 33 






34 


.161458 5.819103 


2 


1 


11 


-.003906 0 56 






35 


.161458 4.252421 


2 


1 


11 


-. 003906 0 57 






36 


.161458 3.469080 


2 


1 


11 


-.003906 0 58 






37 


.161458 2.909551 


2 


1 


11 


-. 003906 059 






38 


.161458 2.461928 


2 


1 


11 


-. 003906 0 60 






39 


,161458 2,238116 


2 


1 


11 


-.003906 0 61 






40 


.161458 2.014305 


2 


1 


11 


-.003906 0 62 






41 


.161458 1.790493 


2 


1 


11 


-.003906 0 63 






42 


.161458 1.342870 


2 


1 


11 


-.003906 0 64 






43 


.161458 0.783341 


2 


1 


11 


-.003906 0 65 






44 


.161458 0 2 1 






11 


-.003906 0 66 
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67 . 147135 5.819103 2 
11 -.0110675 0 89 

68 .147135 4.252421 2 
11 -.0110675 0 90 

69 .147135 3.469080 2 
11 -.0110675 0 91 

70 .147135 2.909551 2 
11 -.0110675 0 92 

71 .147135 2.461928 2 
11 -.0110675 0 93 

72 .147135 2.238116 2 
11 -.0110675 0 94 

73 .147135 2 . 014305 2 
11 -.0110675 0 95 

74 . 147135 1. 790493 2 
11 -.0110675 0 96 

75 . 147135 1. 342870 2 
11 -.0110675 0 97 

76 . 147135 0.783341 2 
11 -.0110675 0 98 

77 .147135 0 2 1 
11 -.0110675 0 99 
00000000000 
ELCON 

I 12 13 2 1 1 1 
1 1 1 1 1 0 10 

II 23 24 13 12 1 1 
1 1 1 1 1 0 20 

21 34 35 24 23 1 1 

I 1 1 1 1 0 30 

31 45 46 35 34 1 1 

II 11 1 0 40 

41 56 57 46 45 1 1 
1 1 1 1 1 0 50 
51 67 68 57 56 1 1 
1 1 1 1 1 0 60 
61 78 79 68 67 1 1 
1 1 1 1 1 0 70 
71 89 90 79 78 1 1 
1 1 1 1 1 0 80 
00000000 



0 0 0 
CTEHP 
1 6 1 
89 99 
0 0 0 
PROP 
1 2 
.1 
4 1 
6 2 

7 2 

8 2 

9 2 

10 2 
PROF 
SYMC 
FORM 
CALC 
PTEMP 
READ, 5 



160 
1 60 
0 0 0 



0 0 0 



.6 70 ( 
0 

10 60 
10 60 
10 60 
10 60 
2 10 60 



1 

1 

1 

1 

1 

1 

1 

1 

1 

1 



0 0 0 0 0 0 



0 0 0 0 
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APPENDIX C 



EXAMPLE TERMINAL SESSION 



The following terminal session was recorded using the 
data file in Appendix B. 



ENTER TERMINAL CODE: 

1 = PLOT- 10 Compatible Terminal (GRAPHICS) 

2 = IBM 3277 DOAL SCREEN (GRAPHICS) 

3 = Any Alpha Numeric Terminal (NO GRAPHICS) 



CP TERMINAL LINESIZE 132 

GLOBAL TXTLIB CMSLIE FORTMO D2 MOD2EEH NONIMSL IMSLSP 
FILEDEF 01 DISK FILE 01 (RECFM VS LRECL 3408 BLKSIZE 34 12)^ 
FILEDEF 08 DISK FILE FT08F001 (RECFM FBA LRECL 132 BLKSIZE 
FILEDEF 04 DISK 

FILEDEF 13 DISK CAL TEST (RECFM FBA LRECL 132 BLKSIZE 132) 
FILEDEF 50 DISK HLP CAL C 

LOAD CALO CAL FRTCHX GROOP1 GROOP2 GROUP3 GROUP4 GR (CLEAR 
EXECUTION BEGINS... 




0.0 SECONDS 



START 

sS**START 



0.003 SECONDS 



PROMPT 

^>»PEOMPT 



0.003 SECONDS 



READ, 4 
4«*READ,4 



0.003 SECONDS 



’^’^HTXFR 



99 ROWS 

1 ROWS 

2 ROWS 
4 ROWS 
1 ROWS 
4 ROWS 

4 ROWS 
7 ROWS 

10 ROWS 

1 ROWS 

2 ROWS 

5 ROWS 
1 ROWS 

9 9 ROWS 
99 ROWS 



1 COLUMNS 
8 COLUMNS 
4 COLUMNS 
1 COLUMNS 
4 COLUMNS 
1 COLUMNS 
4 COLUMNS 
1 COLUMNS 
1 COLUMNS 
99 COLUMNS 
99 COLUMNS 
80 COLUMNS 
99 COLUMNS 
1 COLUMNS 
1 COLUMNS 
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NOMBER OF NODAL POINTS = 99 
NDMBER OF ELEMENTS = 80 
NUMBER OF MATERIAL SETS = 1 
DIMENSION OF COORDINATE SPACE= 2 
DEGREES OF FREEDOM/NODE = 1 
NODES PER ELEMENT (MAXIMUM) = 4 



>S‘>^COORD 






NODE NO 
1 


0. lei 


0.01^ 


2 


0. 166 


0.012 


3 


0. 166 


0.010 


4 


0. 166 


0.008 


5 


0. 167 


0.007 


6 


0. 167 


0.007 


7 


0. 167 


0.006 


8 


0. 167 


0.005 


9 


0. 167 


0.004 


10 


0. 167 


0.002 


11 


0. 167 


0.0 


12 


0.165 


0.017 


13 


0. 165 


0.012 


14 


0, 165 


0.010 


15 


0. 165 


0.008 


16 


0. 165 


0.007 


17 


0. 165 


0.006 


18 


0. 165 


0.006 


19 


0. 165 


0.005 


20 


0. 165 


0.004 


21 


0. 165 


0.002 


22 


0. 165 


0.0 


23 


0. 163 


0.017 


24 


0. 164 


0.012 


25 


0.164 


0.010 


26 


0. 164 


0.008 


27 


0. 164 


0.007 


28 


0. 164 


0.006 


29 


0. 164 


0.006 


30 


0. 164 


0.005 


31 


0. 164 


0.004 


32 


0.164 


0.002 


33 


0. 164 


0.0 


34 


0. 16 1 


0.016 


35 


0. 161 


0.012 


36 


0. 161 


0.010 


37 


0. 161 


0.008 


38 


0. 16 1 


0.007 


39 


0. 161 


0.006 


40 


0. 16 1 


0.006 


41 


0. 16 1 


0.005 


42 


0. 16 1 


0.004 


43 


0. 161 


0.002 


44 


0. 16 1 


0.0 


45 


0. 157 


0.016 


46 


0. 157 


0.012 


47 


0. 157 


0.010 


48 


0. 157 


0.008 


49 


0. 157 


0.007 


50 


0. 157 


0.006 


51 


0. 157 


0.006 


52 


0. 157 


0.005 


53 


0.158 


0.004 


54 


0. 158 


0.002 


55 


0. 158 


0.0 


56 


0. 153 


0.016 


57 


0. 153 


0.011 


58 


0. 153 


0.009 



SECONDS 
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59 

60 
61 
62 

63 

64 

65 

66 

67 

68 

69 

70 

71 

72 

73 

74 

75 

76 

77 

78 

79 

80 
81 
82 

83 

84 

85 

86 

87 

88 

89 

90 

91 

92 

93 

94 

95 

96 

97 

98 

99 



0 . 153 
0 . 154 
0 . 154 
0 . 154 
0 . 154 
0 . 154 
0 . 154 
0 . 154 
0 . 146 
0 . 147 
0 . 147 
0 . 147 
0 . 147 
0 . 147 
0 . 147 
0 . 147 
0 . 147 
0 . 147 
0 . 147 
0 . 135 
0 . 136 
0 . 136 
0 . 136 
0 . 136 
0 . 136 
0 . 136 
0 . 136 
0 . 136 
0 . 136 
0 . 136 
0 , 124 
0 . 125 
0 . 125 
0 . 125 
0 . 125 
0 . 125 
0 . 125 
0 . 125 
0 . 125 
0 . 125 
0 . 125 



'EICON 


f 

1 


2 


3 


4 


5 


1 


12 


13 


2 


1 


1 


2 


13 


14 


3 


2 


1 


3 


14 


15 


4 


3 


1 


4 


15 


16 


5 


4 


1 


5 


16 


17 


6 


5 


1 


6 


17 


18 


7 


6 


1 


7 


18 


19 


8 


7 


1 


8 


19 


20 


9 


8 


1 


9 


20 


2 1 


10 


9 


1 


10 


21 


22 


1 1 


10 


1 


1 1 


23 


24 


13 


12 


1 


1 2 


24 


25 


14 


13 


1 


13 


25 


26 


15 


14 


1 


1 4 


26 


27 


16 


15 


1 


15 


27 


28 


17 


16 


1 


16 


28 


29 


18 


17 


1 


17 


29 


30 


19 


18 


1 


18 


30 


31 


20 


19 


1 


19 


31 


32 


21 


20 


1 


20 


32 


33 


22 


21 


1 


21 


34 


35 


24 


23 


1 


22 


35 


36 


25 


24 


1 


23 


36 


37 


26 


25 


1 



0.008 

0.007 

0.006 

0.005 

0.005 

0.004 

0.002 

0.0 

0.015 

0.011 

0.009 

0.007 

0.006 

0.006 

0.005 

0.005 

0.003 

0.002 

0.0 

0.014 

0.010 

0.008 

0.007 

0.006 

0.005 

0.005 

0.004 

0.003 

0.002 

0.0 

0.013 

0.009 

0.008 

0.006 

0.005 

0.005 

0.004 

0.004 

0.003 

0.002 

0.0 

0 . 083 



SECONDS 



82 



24 


37 


38 


27 


26 


1 


25 


38 


39 


28 


27 


1 


26 


39 


40 


29 


28 


1 


27 


40 


41 


30 


29 


1 


28 


41 


42 


31 


30 


1 


29 


42 


43 


32 


31 


1 


30 


43 


44 


33 


32 


1 


31 


45 


46 


35 


34 


1 


32 


46 


47 


36 


35 


1 


33 


47 


48 


37 


36 


1 


34 


48 


49 


38 


37 


1 


35 


49 


50 


39 


38 


1 


36 


50 


51 


40 


39 


1 


37 


51 


52 


41 


40 


1 


38 


52 


53 


42 


41 


1 


39 


53 


54 


43 


42 


1 


40 


54 


55 


44 


43 


1 


4 1 


56 


57 


46 


45 


1 


42 


57 


58 


47 


46 


1 


43 


58 


59 


48 


47 


1 


44 


59 


60 


49 


48 


1 


45 


60 


61 


50 


49 


1 


46 


61 


62 


51 


50 


1 


47 


62 


63 


52 


51 


1 


48 


63 


64 


53 


52 


1 


49 


64 


65 


54 


53 


1 


50 


65 


66 


55 


54 


1 


51 


67 


68 


57 


56 


1 


52 


68 


69 


58 


57 


1 


53 


69 


70 


59 


58 


1 


54 


70 


71 


60 


59 


1 


55 


71 


72 


61 


60 


1 


56 


72 


73 


62 


61 


1 


57 


73 


74 


63 


62 


1 


58 


74 


75 


64 


63 


1 


59 


75 


76 


65 


64 


1 


60 


76 


77 


66 


65 


1 


6 1 


78 


79 


68 


67 


1 


62 


79 


80 


69 


68 


1 


63 


80 


81 


70 


69 


1 


64 


81 


82 


71 


70 


1 


65 


82 


83 


72 


71 


1 


66 


83 


84 


73 


72 


1 


67 


84 


85 


74 


73 


1 


68 


85 


86 


75 


74 


1 


69 


86 


87 


76 


75 


1 


70 


87 


88 


77 


76 


1 


71 


89 


90 


79 


78 


1 


72 


90 


91 


80 


79 


1 


73 


91 


92 


81 


80 


1 


74 


92 


93 


82 


81 


1 


75 


93 


94 


83 


82 


1 


76 


94 


95 


84 


83 


1 


77 


95 


96 


85 


84 


1 


78 


96 


97 


86 


85 


1 


79 


97 


98 


87 


86 


1 


80 


98 


99 


88 


87 


1 



«^*CTEMP 

NODE 

1 

2 

3 

4 

5 

6 



TEMPERATURE 

160.00 

160,00 

160,00 

160.00 

160.00 

160.00 



0 . 037 



SECONDS 



83 



NODE 


TEMPE RATURE 


89 


60.000 


90 


60.000 


91 


6 0.00 0 


92 


6 0. 000 


93 


60.000 


94 


60.000 


95 


6 0.00 0 


96 


60.000 


97 


6 0.000 


98 


6 0.00 0 


99 


60.000 


^SjspHOP 



0.007 SECONDS 



MATERIAL SET 1 FOR ELEMENT TYPE 2 
DEGREE OF FREEDOM ASSIGNMENTS 

HEAT CONDUCTION ELEMENT 

CONDUCTIVITY KX = .1000000 

SPECIFIC HEAT .6000000 

DENSITY 70, 00000 

HEAT GENER/UNIT VOL .0 

4 GAUSS PT S/DIR 

5 LINES WITH SPECIFIED BOUNDARY 
PLANE ANALYSIS 

5 ROWS 3 COLUMNS 

3 ROWS 9 COLUMNS 

5 ROWS 2 COLUMNS 

LINE B.C. 

ELEM B.C. Lip PROgEgJY^ VALUE 

7 2 2 lOl 00000 

8 2 2 10 . 00000 

9 2 2 10. 00000 

10 2 2 10. 00000 

0.027 

si«*!«PROF 

82 ROWS 1 COLUMNS 

0.017 

=3»S!S YMC 

945 ROWS 1 COLUMNS 

0.406 

!»$F0RM 

FLUX CONVERGENCE TEST 

RNMAX = 60.927 RN 

0.716 

$«sCALC 

ENERGY (DR^»A«^DR) = 0. 2 2491 7705 2D + 0 5 

L 0.030 

‘S’S'ETEMP 



LOCAL GLOBAL 

NUMBER NUMBER 

1 1 



KY = . 1000 



CONDITIONS 



60.00000 
60. 00000 
60.00000 
60.00000 
SECONDS 



SECONDS 



SECONDS 



= 60.927 

SECONDS 



SECONDS 



NODAL TEMPERATURES 



TIME 



NODE 

1 

2 

3 

4 

5 

6 
7 



TEMP 

0. 16000D+03 
0. 16000D+03 
0. 16000D+03 
0. 16000D+03 
0. 16000D+03 
0. 16000D+03 
0. 13970D+03 



0.0 
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8 


0. 1299 5D+03 


9 


0. 11974D+03 


10 


0. 1136 2D+03 


1 1 


0. 11094D+03 


12 


0. 15581D+03 


13 


0.15526D+03 


14 


0.1541 9D+03 


15 


0.15233D+03 


16 


0. 14821D+03 


17 


0. 14317D+03 


18 


0. 13992D+03 


19 


0. 13434D+03 


20 


0. 1256 8D+03 


21 


0. 11957D+03 


22 


0. 11 680D+03 


23 


0. 15 163D+03 


24 


0. 15064D+03 


25 


0. 14880D+03 


26 


0. 14597D+03 


27 


0. 14199D+03 


28 


0. 13 96 2D+03 


29 


0. 13664D+03 


30 


0. 13384D+03 


31 


0. 12838D + 03 


32 


0. 1234 9D+03 


33 


0.12 104D+03 


34 


0. 1435 5D+03 


35 


0. 1421 4D+03 


36 


0. 13999D+03 


37 


0. 13763D+03 


38 


0. 13 53 0D+03 


39 


0. 13400D+03 


40 


0. 13274D+03 


41 


0. 13142D+03 


42 


0. 12897D+03 


43 


0. 12656D+03 


44 


0.12520D+03 


45 


0. 13257D+03 


46 


0. 13148D+03 


47 


0.13016D+03 


48 


0. 12898D+03 


49 


0. 12795D+03 


50 


0. 1274 4D+03 



NODAL TEHPEEATORES TIME 0.0 

TEMP 

0.12693D+03 
0. 12644D+03 
0. 12556D+03 
0. 12 47 3D+03 
0. 12427D+03 
0. 12 30 0D+03 
0. 12244D+03 
0. 12181D+03 
0, 12127D+03 
0. 12084D+03 
0. 1206 2D+03 
0. 1204 ID+03 
0. 1202 1D+03 
0. 1198 7D+03 
0. 11955D+03 
0. 1193 7D+03 
0. 10875D+03 
0. 1086 2D+03 
0. 1084 8D+03 
0. 10 83 6D+03 



NODE 

51 

52 

53 

54 

55 

56 

57 

58 

59 

60 
61 
62 

63 

64 

65 

66 

67 

68 

69 

70 
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71 


0. 10827D+03 


72 


0. 1082 3D+03 


73 


0. 1081 9D+03 


74 


0. 1081 5D+03 


75 


0. 10808D+03 


76 


0.10801D+03 


77 


0. 10798D+03 


78 


0. 8516 2D+02 


79 


0. 8516 3D+02 


80 


0.85163D+02 


81 


0.85163D+02 


82 


0.85 163D+02 


83 


0.85163D+02 


84 


0.85163D+02 


85 


0.85163D+02 


86 


0.85163D+02 


87 


0.8516 3D+02 


88 


0.85163D+02 


89 


0.60 00 0D+02 


90 


0.60000D+02 


91 


0.60000D+02 


92 


0.60000D+02 


93 


0. 60 00 0D+ 02 


94 


0.60000D+02 


95 


0.60000D+02 


96 


0.60 00 0D+02 


97 


0.60 00 0D+02 


98 


0.60000D+02 


99 


0. 60000D+02 


5 



Stop 

JStjtsTOP 

R; T= 1.80/2.67 12:36: 15 



0.043 SECONDS 
0.003 SECONDS 



f 
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